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Purposes Transforming growth factor pi (TGF-pi) inhibits the growth of bladder cancer 
cells and this effect is prominent and constant in 253 J bladder cancer cells. We per- 
formed a microarray analysis to search for genes that were altered after TGF-pi treat- 
ment to understand the growth inhibitory action of TGF-pi. 

Materials and Methods: 253J bladder cancer cells were exposed to TGF-pl and total 
RNA was extracted at 6, 24, and 48 hours after exposure. The RNA was hybridized onto 
a human 22K oligonucleotide microarray and the data were analyzed by using 
GeneSpring 7.1. 

Results: In the microarray analysis, a total of 1,974 genes showing changes of more than 
2.0 fold were selected. The selected genes were further subdivided into five highly cohe- 
sive clusters with high probability according to the time-dependent expression pattern. 
A total of 310 genes showing changes of more than 2.0 fold in repeated arrays were iden- 
tified by use of simple t-tests. Of these genes, those having a known function were listed 
according to clusters. Microarray analysis showed increased expression of molecules 
known to be related to Smad-dependent signal transduction, such as SARA and Smad4, 
and also those known to be related to the mitogen-activated protein kinase (MAPK) 
pathway, such as MAPKK1 and MAPKK4. 

Conclusions: A list of genes showing significantly altered expression profiles after 
TGF-pi treatment was made according to five highly cohesive clusters. The data sug- 
gest that the growth inhibitory effect of TGF-pi in bladder cancer may occur through 
the Smad-dependent pathway, possibly via activation of the extracellular signal-re- 
lated kinase 1 and Jun amino-terminal kinases Mitogen-activated protein kinase 
pathway. 
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INTRODUCTION 

Transforming growth factor p (TGF-p) is a member of a fam- 
ily of dimeric polypeptide growth factors that includes bone 
morphogenic proteins and activins [1]. Every cell in the 
body, including epithelial, endothelial, hematopoietic, 
neuronal, and connective-tissue cells, produces TGF-p and 
has receptors for it [2]. Inhibition of cell proliferation is cen- 



tral to the TGF-p response in the epithelial lineage and es- 
cape from this response is a hallmark of many cancer cells 
[3]. 

In contrast with these opinions, we found that most blad- 
der cancer cell lines are sensitive to the growth inhibitory 
action of TGF-pi [4] . If human bladder cancers exhibit sen- 
sitivity to the growth inhibitory action of TGF-pi, TGF-pi 
may be a strong candidate molecule for treating this hor- 
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rible disease. 

In a previous study, we found that TGF-pi inhibits the 
cellular growth of several bladder cancer cell lines. We thus 
assumed that bladder cancer cells are sensitive to the 
growth inhibitory action of TGF-pi. The aim of this study 
was to investigate how TGF-p 1 inhibits the cellular growth 
of bladder cancer cells. For this, we investigated altered 
gene expression profiles after TGF-pi treatment in 253J 
bladder cancer cells. We compared the altered gene ex- 
pression profiles obtained with the currently suggested 
signal transduction pathway of TGF-p to infer the mecha- 
nism of the growth inhibition of 253 J bladder cancer cells 
by TGF-pi. 

MATERIALS AND METHODS 

1. Cells and culture conditions 

The human bladder cancer cell line 253J was obtained from 
the Korea Cell Line Bank (Seoul National University, 
Seoul, Korea). The cells were maintained in Dulbecco's 
modified Eagle medium containing 10% fetal bovine se- 
rum, 100 units of penicillin/mL, and 100 ug of streptomy- 
cin/mL. TGF-pi was purchased from Sigma Chemical Co. 
(St. Louis, MO, USA). 

2. RNA extraction 

Total RNA was isolated from cells that were grown to ap- 
proximately 60% confluence in 250-mL culture flasks 
(Sigma Chemical Co.) by use of TRI reagent (Gibco 
BRL/Life Technologies, Grand Island, NY, USA). The total 
RNA was phenol/chloroform-extracted, ethanol precipita- 
ted, and cleaned with RNeasy cleanup system columns 
(Qiagen, Valencia, CA, USA). The quantity and quality 
were determined by optical density measurements at 260 
and 280 nm. 

3. Microarray analysis 

The human 22K oligonucleotide chip (Illumina Oligonu- 
cleotide Library, San Diego, CA, USA) was used in this 
study. Each 10 ug of total RNA was reverse transcribed in 
the presence of Cy3- or Cy5-dUTP (NEN Life Sciences, 
Boston, MA, USA) at 42°C for 2 hours. Control RNA was 
labeled with fluorescent Cy3-dUTP and test condition RNA 
was labeled with fluorescent Cy5-dUTP. Both the Cy3- and 
Cy5-labeled cDNA were purified by using the polymerase 
chain reaction (PCR) purification kit (Qiagen, Valencia, 
CA, USA) as recommended by the manufacturer. The puri- 
fied cDNA was resuspended in 100 uL of hybridization sol- 
ution containing 5x saline-sodium citrate (SSC), 0.1% so- 
dium dodecyl sulfate (SDS), 30% formamide, 20 ug of 
Human Cot-1 DNA, 20 ug of poly A RNA, and 20 ug of yeast 
tRNA (Invitrogen, Grand Island, NY, USA). The hybrid- 
ization mixtures were heated at 100°C for 2 to 3 minutes 
and were directly pipetted onto microarrays. The arrays 
were allowed to hybridize at 42°C for 12 to 16 hours in the 
humidified hybridization chamber (GenomicTree Inc., 
Daejeon, Korea). The hybridized microarrays were washed 



with 2x SSC/0.1% SDS for 5 minutes, O.lx SSC/0.1% SDS 
for 10 minutes, and O.lx SSC for 2 minutes two times. The 
washed microarrays were immediately dried by using the 
microarray centrifuge (GenomicTree Inc., Daejeon, Korea). 

4. Microarray data acquisition 

The hybridization images were analyzed by using GenePix 
Pro 4.0 (Axon Instruments, Union City, CA, USA). The 
average fluorescence intensity for each spot was calculated 
and local background was subtracted. All data normal- 
ization and selection of fold-changed genes were performed 
by using GeneSpring 7.1 (Silicon Genetics, Redwood City, 
CA, USA). The reliable genes were filtered with a cutoff val- 
ue based on the two-component error model after in- 
tensity-dependent normalization (LOWESS) [5,6]. The 
averages of the normalized ratios were calculated by divid- 
ing the average of normalized signal channel intensity by 
the average of normalized control channel intensity. The 
analysis of variance test (parametric) and single t-test were 
performed at p-values of < 0.01 or < 0.05 to find genes that 
were differentially expressed across conditions. Unsuper- 
vised hierarchical clustering was performed by similarity 
measurements based on Pearson correlations around zero. 
Functional annotation of genes was performed according 
to the Gene Ontology Consortium (http://www.geneontolo- 
gy.org/index.shtml) by use of GeneSpring 7.1. 

5. Expression of selected genes by RT-PCR 

Total RNA from the Trizol isolate was treated with 
RNase-free DNase I. After removal of the DNase 1, 1 ug of 
total RNA was reverse-transcribed by using random hex- 
adeoxynucleotide primers. The genes of interest, Samd4, 
Gadd45p, TGF-pi, and the housekeeping gene p-actin, 
were analyzed by use of a Quantum RNA reverse tran- 
scription PCR (RT-PCR) kit according to the manu- 
facturer's protocol (Ambion, Austin, TX, USA). The follow- 
ing primers were used to amplify Smad4: forward, 
5'-CCCAGGATCAGTAGGTGGAA-3', and reverse, 5'-CC- 
ATGC CTGAC AAGTTCTGA (451-bp); Gadd45p. forward, 
5'-GGAAGGTTTTGGGCTCTCTG-3', and reverse, 5'-GT- 
GTGAGGGTTCGTGACCAG-3' (471-bp); Smad3: for- 
ward, 5'-GGGCTCCCTCATGTCATCTA-3', and reverse, 
5'-GGCTCGCAGTAGGTAACTGG-3' (441-bp); TGF-ft: 
forward, 5'-GCCCTGAGAGACCTGCTGAA-3', and re- 
verse, 5'-TCGCCTTCCCGTTGATAGTG-3' (147-bp). Beta- 
actin primers (Ambion) were used as an internal standard 
(294 bp). PCR was performed as 94°C for 30 seconds, 56°C 
for 30 seconds, and 72°C for 30 seconds for 31 cycles, fol- 
lowed by a final elongation for 7 minutes. PCR products 
were electrophoresed on a 1.5% agarose gel and were vi- 
sualized by ethidium bromide staining. 

RESULTS 

1. Gene expression patterns after TGF-p treatment in the 
253 J cell line 

In the microarray analysis, the gene expression patterns 
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FIG. 1. Hierarchical clustering of the gene expression profiles of 7,714 reliable genes (A) and 1,974 genes showing changes ^2.0 fold 
in at least 1 array (B). Gene ontology map (C) and classification of 1,974 reliable genes according to GeneSpring 7.1 (D). 



were observed at 6, 24, and 48 hours after TGF-p treatment. 
The genes were divided into TGF-p 1-induced, TGF-pi-re- 
pressed, and other genes from 6 hours after the TGF-pi 
treatment. The expression pattern was persistent and in- 
tensified for up to 48 hours after the TGF-pi treatment. We 
filtered the analysis for reliable genes that had a sum of me- 
dian > 1,500 in 3 arrays (sum of median=Cy3 net sig- 
nal+Cy5 net signal). This resulted in 7,714 genes. 

Hierarchical cluster analysis was used to profile the gene 
expression patterns (Fig. 1A). Those genes showing 
changes of more than 2.0 fold in at least 1 array were 
selected. A total of 1,974 genes were found (Fig. IB). These 
genes were classified by their known function (Fig. 1C, D). 
Because the expression difference between TGF-p treat- 
ment and the control was most remarkable at 48 hours after 
TGF-p treatment, we repeated the microarray at 48 hours. 
In this analysis, 7,992 reproducible genes were filtered, 
and a total of 310 genes showing similar expression in the 
two microarrays were selected by use of simple t-tests (Fig. 



2). We then categorized these genes according to their 
known function (Table 1). 

2. Cluster analysis 

The genes showing changes of more than 2.0 fold in at least 
1 array totaled 1,974 genes. Among them, a large number 
of genes were further subdivided into 5 highly cohesive 
clusters with high probability according to the time-de- 
pendent expression pattern (Fig. 3). After excluding genes 
whose name or function was not yet reported (i.e., hypo- 
thetical proteins), 270 genes were included in the cluster 
analysis. Other significant gene lists of 310 genes were ob- 
tained by using the two microarrays and simple t-tests (Fig. 
2). 

We selected 310 genes and performed a cluster analysis 
and made a gene list according to five different clusters. The 
gene profiles of each cluster are shown in Tables 2-6. 
Cluster 1 genes showed increased expression over the 
whole time period but were more up-regulated after 24 
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Hierarchical clustering 
of reliable genes 
(7,992 genes ) 
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Hierarchical clustering 
of statistically significant 
genes (T test, p<0.01) 

(310 genes ) 



FlG. 2. Hierarchical clustering of the gene expression profiles of 
7,992 reliable genes (A) and 310 statistically significant genes 
(B). 



Table 1. Classification of 310 statistically significant genes as 
filtered by a repeat microarray 



Group 


Description 


No. of genes 


1 


Apoptosis regulator 


8 


2 


Cancer genes 


3 


3 


Cell cycle regulator 


7 


4 


Chaperone 


8 


5 


Enzyme 


72 


6 


Immunity protein 


3 


7 


Microtubular dynamics 


1 


8 


Motor 


1 


9 


Nucleic acid binding 


35 


10 


Signal transducer 


10 


11 


Storage 


0 


12 


Structural protein 


32 


13 


Transport 


34 


14 


Unclassified 


160 



hours sequentially. Cluster 2 genes showed decreased ex- 
pression over the whole time period. Cluster 3 genes 
showed meaningful changes at 6 or 24 hours but became 
up-regulated at 48 hours. Cluster 4 genes did not show 
meaningful changes at 6 or 24 hours but were down-regu- 
lated at 48 hours. Cluster 5 genes showed markedly in- 
creased expression over the whole time period but were no- 
ticeably up-regulated at 48 hours after TGF-p 1 treatment. 

3. RT-PCR analysis 

RT-PCR was used to evaluate a subset of genes identified 
by microarray analysis as undergoing significant changes 
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Fig. 3. Cluster analysis of the 310 statistically significant genes 
according to the time response. These genes were classified into 
5 clusters by GeneMaths. Cluster 1 (red); showing increased 
expression over the whole time periods but more up-regulated 
after 24 hours sequentially, Cluster 2 (green); showing 
decreased expression over the whole time periods, Cluster 3 
(purple); without showing meaningful changes at 6 or 24 hours 
but becoming up-regulated at 48 hours, Cluster 4 (blue); without 
showing meaningful changes at 6 or 24 hours but becoming 
down-regulated at 48 hours, Cluster 5 (yellow); showing 
markedly increased expression over the whole time periods but 
becoming noticeably up-regulated only at 48 hours after 
transforming growth factor |31 treatment. 



in expression. We selected three interesting genes for 
RT-PCR confirmation: Samd4, Gadd45p, and TGF-pi. 
Specific primers were designed for the selected genes and 
mRNA expression was checked by RT-PCR with the use of 
P-actin as an internal standard. The expression of these 
genes in RT-PCR was comparable to the results obtained 
by microarray (Fig. 4). 



Korean J Urol 2014;55:542-550 



546 



Lee et al 



TABLE 2. Cluster 1: a total of 27 genes showed a steady increase in expression over the whole time range 



Gene name 


Gene bank 
accession no. 


Protein description ~ 


6h 


Fold change 
24 h 


48 h 


RAD54B 


NM. 


.006550 


RAD54B 


1.5 


1.5 


3.6 


CBX1 


NM. 


.006807 


Chromobox homolog 1 


1.5 


1.3 


2.4 


DERI 


NM. 


.018630 


Derlin-1 


1.4 


1.4 


2.5 


RBBP4 


NM. 


.005610 


Retinoblastoma binding protein 4 


1.4 


1.4 


2.2 


ITGB3BP 


NM. 


.014288 


Integrm beta 3 binding Protein 


1.3 


1.2 


2.2 


RIF1 


NM. 


.018372 


Receptor-interacting factor 1 


1.2 


1.4 


2.2 


CN0T7 


NM. 


.013354 


CCR4-NOT transcription complex, subunit 7 


1.5 


1.4 


2.4 


PSMA2 


NM. 


.002787 


Proteasome subunit, alpha type, 2 


1.4 


1.4 


2.2 


LEPRE1 


NM. 


.022356 


Leucine proline-enriched Proteoglycan 


1.1 


1.3 


2.3 


GTF2H1 


NM. 


.005316 


General transcription factor IIH 


1.5 


1.3 


2.2 


KIAA0934 


XM_ 


034536 


catalytic activity 


1.1 


1.1 


2.1 


SEMA3C 


NM. 


.006379 


r~i l •• 111*1 •iii , i * ii 

bema domain, immunoglobulin domain, short basic domain, secreted, 


1.2 


1.4 


3.0 








(semaphorin) 3C 








sun 


NM. 


005801 


Putative translation initiation factor 


1.2 


1.3 


2.7 


C130RF7 


NM. 


.024546 


Chromosome 13 open reading frame 7 


1.2 


1.1 


2.4 


ETFA 


NM^ 


.000126 


Electron-transfer-flavoprotein, alpha polypeptide 


1.3 


1.2 


2.6 


C120RF14 


NM. 


.021238 


Chromosome 12 open reading frame 14 


1.5 


1.6 


3.0 


EPHB4 


NM. 


004444 


EphB4 


1.1 


1.2 


2.2 


LOC54499 


NM. 


.019026 


Putative membrane protein 


1.4 


1.4 


2.5 


PITPNB 


NM. 


.012399 


Phosphatidylinositol transfer protein, beta 


1.3 


1.1 


2.0 


KIAA1078 


XM_ 


036589 


Homo sapiens KIAA1078 Protein 


1.7 


1.6 


3.5 


COPA 


NM. 


.004371 


Coatomer protein complex, subunit alpha 


1.2 


1.7 


2.4 


CD44 


NM. 


.000610 


CD44 antigen, homing function and Indian blood group system 


1.4 


1.9 


3.9 


HECTD1 


XM_ 


030175 


HECT domain containing 1 


1.1 


1.0 


2.0 


C210RF7 


NM. 


.020152 


Chromosome 21 open reading frame 7 


1.2 


1.4 


2.3 


UFM1 


NM. 


.016617 


Ubiquitin-fold modifier 1 


1.1 


1.3 


2.1 


TULP4 


IN 1V1 


0909/1 K 
UZUZ40 


Tubby like protein 4 


1.7 


1.6 


3.4 


SH3KBP1 


NM. 


031892 


SH3-domain kinase binding protein 1 


1.2 


1.3 


2.0 


MPZL1 


NM. 


_024569 


Myelin protein zero-like 1 


1.0 


1.2 


2.5 


KIF21A 


XM_ 


.040211 


Kinesin familv memher 21 A (KTF21 A) mRNA 

_L VlllvT o 111 IcLlllll V 111CT111 KJKZl. £j _LiTi \ _L VI. i 111/ ^ 111XL1. > H> 


1.3 


1.1 


2.2 


APLP2 


NM. 


_001642 


Amyloid beta (A4) precursor-like protein 2 


1.2 


1.4 


2.8 


SMBP 


NM. 


020123 


SM- 11044 binding protein 


1.1 


1.2 


2.2 


PDCD4 


NM. 


_145341 


Programmed cell death 4, neoplastic transformation inhibitor 


1.1 


1.5 


2.2 


NUDT4 


NM. 


_019094 


Nucleoside diphosphate linked moiety X-type motif 4 


1.1 


1.1 


2.4 


FAT 


NM. 


.005245 


FAT tumor suppressor homolog 1 (Drosophila) 


1.4 


1.3 


2.2 


LD0C1L 


NM. 


_032287 


Leucine zipper, down-regulated in cancer 1-like 


1.2 


1.5 


2.1 


RAB10 


NM^ 


_016131 


Member RAS oncogene family 


1.1 


1.2 


2.0 


HNRPK 


NM. 


_002140 


Heterogeneous nuclear ribonucleoprotein K 


1.4 


1.3 


2.0 



4. Assumption of the growth inhibitory pathway of TGF-pi 
on 253 J cells 

1) Smad-dependent pathway 

Among 1,974 genes showing changes of more than 2.0 fold 
in at least one array, four representative Smad-dependent 
pathway-related genes were detected. These genes in- 
cluded TGIF, Smad4, Smad5, and SARA1. Microarray 
analysis showed increased expression of SARA and the 
common Smad, Smad4 (Fig. 5). This result indicated that 
the growth inhibitory signal transduction of TGF-|31 in 
bladder cancer might be mediated by the Smad-dependent 
pathway. 



2) Smad-independent pathway 

Among the 1,974 genes, seven Smad-independent mi- 
togen-activated protein kinase (MAPK) pathway path- 
way-related genes were detected (Fig. 6). Of these, 
MAPKK1 and MAPKK4 expression were up-regulated. 
MAPKK1 is known to activate extracellular signal-related 
kinase 1 (ERK 1), and MAPKK4 and MAPKKK7 are known 
to activate Jun amino-terminal kinases (JNK). This result 
indicated that the growth inhibitory signal transduction of 
TGF-pl in bladder cancer might be mediated by the 
Smad-independent MAPK pathway, especially the ERK 1 
and JNK pathways. 
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Gene name 


Gene bank 
accession no. 


Protein description 


C v. 
d n 


Fold change 

9 A Vi 

Z^i 11 


*±o 11 


C4BPA 


NM. 


000715 


Complement component 4 binding protein, alpha 


0.2 


0.4 


0.5 


GAN 


NM. 


.022041 


Giant axonal neuropathy (gigaxonin) 


0.9 


0.6 


0.4 


CASP8AP2 


NM. 


012115 


CASP8 associated protein 2 


0.2 


0.2 


0.3 


ZNF167 


NM. 


.018651 


Zinc finger protein 167 


0.1 


0.2 


0.2 


NLGN4Y 


NM. 


.014893 


Neuroligin 4, Y-linked 


0.4 


0.3 


0.2 


TUBGCP6 


NM. 


.020461 


Tubulin, gamma complex associated protein 6 


0.6 


0.6 


0.4 


SLMAP 


NM. 


.007159 


Sarcolemma associated protein 


0.3 


0.2 


0.2 


RGPR 


NM. 


.033127 


Regucalcin gene promotor region related protein 


0.7 


0.6 


0.5 


HOM-TES-103 


NM. 


.080731 


HOM-TES-103 tumor antigen-like 


0.6 


0.5 


0.6 


PRKCDBP 


NM. 


.145040 


Protein kinase C, delta binding Protein 


0.5 


0.7 


0.7 


SEC23A 


NM. 


.006364 


Sec23 homolog A Saccharomyces cerevisiae 


0.5 


0.6 


0.3 


ZNF43 


NM. 


.003423 


Zinc finger protein 43 (HTF6) 


0.4 


0.7 


0.5 


SLC7A4 


NM. 


004173 


Solutecarrier family 7, member 4 


0.6 


0.8 


0.4 


RNASE6 


NM. 


.005615 


Ribonuclease, RNase A family, k6 


1.0 


0.5 


0.6 


RUNX3 


NM. 


.004350 


Runt-related transcription factor 3 


0.6 


0.7 


0.4 


COL9A1 


NM. 


.078485 


Collagen, type IX, alpha 1 


0.5 


0.7 


0.3 


CLMN 


NM. 


.024734 


Calmin (calponin-like, transmembrane) 


0.7 


0.7 


0.5 


ARHGEF1 


NM. 


.004706 


Rho guanine nucleotide exchange factor (GEF) 1 


0.6 


0.6 


0.4 


HYPE 


NM. 


.007076 


Huntingtin interacting protein E 


0.8 


0.5 


0.4 


ACOX1 


NM. 


.007292 


Acyl-Coenzyme A oxidase 1, palmitoyl 


0.1 


0.5 


0.2 


MMP25 


NM^ 


.004142 


Matrix metalloproteinase-like 1 


0.7 


0.5 


0.4 


REC8L1 


NM. 


.005132 


REC8-like 1 (yeast) 


0.7 


0.9 


0.3 


MAL 


NM. 


.002371 


Mai, T-cell differentiation Protein 


0.8 


0.7 


0.3 


SE57-1 


NM. 


.025214 


CTCL tumor antigen se57-l 


0.3 


0.5 


0.4 


CDKN2A 


NM. 


.058196 


Cyclin-dependent kinase inhibitor 2A 


0.7 


0.7 


0.5 


C220RF1 


NM. 


.001585 


Ch'some 22 open reading frame 1 


0.7 


0.8 


0.3 


SLC24A6 


NM. 


.024959 


Solute carrier family 24, member 6 


0.5 


0.6 


0.3 


SERPINA6 


NM. 


.001756 


Serine proteinase inhibitor, clade A, member 6 


0.5 


0.5 


0.1 


F0LR1 


NM. 


.016730 


Folate receptor 1 (adult) 


0.7 


0.7 


0.4 


CRN7 


NM. 


.024535 


Coronin 7 


0.7 


0.6 


0.4 


EVER1 


NM. 


.007267 


Epidermodysplasia verruciformis 1 


0.7 


0.7 


0.4 


C90RF58 


NM. 


.031426 


Chromosome 9 open reading frame 58 


0.4 


0.4 


0.1 


MGLL 


NM. 


.007283 


Monoglyceride lipase 


0.8 


0.5 


0.4 


LOC81558 


NM. 


.030802 


C/EBP-induced protein 


0.5 


0.3 


0.5 


TAGLN 


NM. 


.003186 


Transgelin 


0.3 


0.3 


0.3 


KCNJ15 


NM. 


.002243 


Potassium inwardly-rectifying channel, subfamily J, member 15 


0.5 


0.6 


0.7 


RAD52 


NM. 


.134422 


RAD52 homolog (S. cerevisiae) 


0.5 


0.4 


0.4 


TPPP 


NM. 


.007030 


Brain-specific protein p25 alpha 


0.4 


0.4 


0.4 


AKIP 


NM. 


017900 


Aurora-A kinase interacting Protein 


0.7 


0.5 


0.4 


KLF1 


NM. 


.006563 


Rruppel-like factor 1 (erythroid) 


0.1 


0.5 


0.1 


SLC30A5 


NM. 


.024055 


Solute carrier family 30 (zinc transporter), member 5 


0.6 


0.8 


0.5 


GNG13 


NM. 


.016541 


Guanine nucleotide binding protein (G protein), gamma 13 


0.6 


0.5 


0.4 


HTR3B 


NM. 


.006028 


5-Hydroxytryptamine (serotonin) receptor 3B 


0.8 


0.3 


0.7 



TABLE 4. Cluster 3: a total of 3 genes showed no meaningful changes at 6 or 24 hours but significant up-regulation at 48 hours 



Gene name 


Gene bank accession no. 


Protein description 




Fold change 




6h 


24 h 


48 h 


GOLPH4 


NM _0 14498 


Golgi phosphoprotein 4 


0.9 


0.9 


2.0 


RPIA 


NM_144563 


Ribose 5-phosphate isomerase A 


1.2 


0.8 


2.0 


PLK2 


NM_006622 


Polo-like kinase 2 (drosophila) 


1.2 


1.0 


2.3 
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TABLE 5. Cluster 4: a total of 12 genes showed no meaningful changes at 6 or 24 hours but became down-regulated at 48 hours 



Gene name 


Gene bank 


Protein description 


Fold chang 


;e 


accession no. 


6h 


24 h 


A Q V, 
4o 11 


PCSK6 


NM. 


138319 


Proprotein convertase subtilisin/kexin type 6 


1.5 


0.7 


0.4 


EGFR 


NM_ 


005228 


Epidermal growth factor receptor 


0.8 


0.9 


0.3 


HCRTR1 


NM_ 


001525 


Hypocretin receptor 1 


1.2 


0.9 


0.5 


PIK3R2 


NM_ 


005027 


Phosphoinositide-3-kinase, regu latory subunit, polypeptide 2 


1.0 


1.2 


0.5 


MMAB 


NM_ 


052845 


Methylmalonic aciduria type B 


0.9 


1.2 


0.5 


CGI-14 


NM_ 


015944 


CGI-14 protein 


1.2 


1.2 


0.5 


ABCF1 


NM_ 


001090 


ATP-binding cassette, subfamily F (GCN20), member 1 


1.0 


1.1 


0.5 


PRSS11 


NM_ 


002775 


Protease, serine, 11 (IGF binding) 


1.0 


0.8 


0.4 


8D6A 


NM_ 


016579 


8D6 antigen 


1.2 


1.3 


0.5 


MAGMAS 


NM_ 


016069 


Mitochondria-associated protein involved in GMCSF signal transduction 


0.9 


1.0 


0.4 


CLN3 


NM_ 


000086 


Ceroid-lipofuscinosis, neuronal3 


0.9 


1.0 


0.5 


AKR1B1 


NM_ 


001628 


Aldo-keto reductase family 1, member Bl 


1.0 


0.8 


0.5 



TABLE 6. Cluster 5: a total of 21 genes reacted with a noticeable increase in expression only at 48 hours 



Gene bank _ i . . t . Fold change 

Gene name . Protein description 

accession no. 6 h 24 h 48 h 



PREI3 


NM_ 


015387 


Preimplantation protein 3 


1.9 


1.9 


2.4 


HEY1 


NM_ 


012258 


Hairy/enhancer-of-split related with YRPW motif 1 


1.4 


2.7 


3.0 


ZDHHC14 


NM_ 


024630 


Zinc finger, DHHC domain containing 14 


1.9 


1.4 


2.0 


GC20 


NM_ 


005875 


Translation factor suil homolog 


1.3 


1.7 


2.1 


GADD45B 


NM_ 


015675 


Growth arrest and DNA-damage-inducible, beta 


2.8 


3.4 


6.8 


DGUOK 


NM_ 


080917 


Deoxyguanosine kinase 


1.3 


1.6 


2.1 


RBPSUH 


NM_ 


015874 


Recombining binding protein suppressor of hairless(Drosophila) 


1.4 


1.5 


2.1 


DDEF2 


NM_ 


003887 


Development and differentiation enhancing factor 2 


1.8 


1.8 


3.2 


CHPPR 


NM_ 


014637 


Likely ortholog of chicken chondrocyte protein with a poly-proline region 


1.4 


1.7 


2.2 


ASE-1 


NM_ 


012099 


CD3-epsilon-associated protein; antisense to ERCC-1 


1.8 


1.7 


2.4 


ITGA6 


NM_ 


000210 


Integrin, alpha 6 


1.4 


1.5 


2.1 


ANGPTL4 


NM. 


_139314 


Angiopoietin-like 4 


5.1 


7.8 


6.8 


NUP153 


NM_ 


005124 


Nucleoporin 153kDa 


1.5 


1.5 


2.3 


SYPL 


NM_ 


006754 


Synaptophysin-like protein 


1.3 


1.6 


2.1 


NDRG1 


NM_ 


006096 


N-myc downstream regulated gene 1 


1.3 


1.7 


2.2 


ADA 


NM_ 


000022 


Adenosine deaminase 


2.0 


3.0 


2.2 


PDLIM7 


NM_ 


005451 


PDZ and LIM domain 7 


1.8 


3.0 


2.9 


SCD 


NM_ 


005063 


Stearoyl-CoA desaturase (delta-9-desaturase) 


1.5 


2.7 


4.6 


SLC12A1 


NM_ 


000338 


Solute carrier family 12, member 1 


2.7 


3.4 


6.3 


HMOX1 


NM_ 


002133 


Heme oxygenase (decycling) 1 


2.2 


2.2 


3.2 


SLC20A1 


NM_ 


005415 


Solute carrier family 20, member 1 


1.7 


2.7 


4.4 



DISCUSSION 

Microarray analysis was performed to search for genes that 
were altered after the TGF-pi treatment to study the mech- 
anisms of the growth inhibitory action of TGF-pi. 253J cells 
were selected because this cell line showed constant growth 
inhibition in the cell viability assays [4]. The 253J cells 
were treated with TGF-pi and were harvested at 6, 24, and 
48 hours after treatment. In the microarray analysis, the 
genes were divided into TGF-pi-induced, TGF-pi-re- 
pressed, and other genes from 6 hours after the TGF-pi 
treatment. A total of 7,714 altered genes were filtered after 
the TGF-pi treatment. Among these genes, the genes 



showing changes of more than 2.0 fold in at least one array 
were selected. A total of 1,974 genes were found. These 
genes were classified along with their known function (Fig. 
1). 

The microarray analysis was repeated by using 48-hour 
RNAs, because the expression pattern was persistent and 
intensified for up to 48 hours after the TGF-pi treatment. 
We filtered 7,992 reproducible genes through the repeat 
microarray and found a total of 310 genes showing similar 
expression in the two microarrays by use of simple t-tests 
(Fig. 2A, B). We categorized these genes according to their 
known function. Seventy-two genes were categorized as 
enzymes. Thirty-five genes were categorized as genes asso- 
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Fig. 5. The expression patterns of the Smad-dependent pathway 
related genes. Most signal transducers showed elevated 
expression. 
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Fig. 6. The expression patterns of the MAP kinsae pathway 
related genes. 



dated with nucleic acid binding, and another 34 genes were 
categorized as genes associated with transport. Thirty-two 
genes were classified in the structural protein category, 8 
genes were associated with apoptosis, 3 genes were cancer 
genes, and 7 genes were associated with the cell cycle (Table 
1). 

We performed a cluster analysis to estimate the function 
of selected genes to gain insight into the mechanism of the 
TGF-pi effect. Searching for meaningful information pat- 
terns in gene expression data is not trivial. An initial step 
was to cluster or group genes with similar changes in ex- 
pression [7]. Moreover, we traced groups of genes showing 
similar expression patterns according to the time of ex- 
pression and subdivided these into 5 highly cohesive clus- 
ters (Fig. 3). The cellular reaction process related to ex- 



ternal stimuli or an altered internal milieu is generally se- 
quential, and a variety of genes or proteins respond 
simultaneously. Thus, the assumption that genes showing 
a similar expression pattern according to reaction time will 
have a similar function may be true. According to this as- 
sumption, these selected 1,974 genes were further sub- 
divided into 5 highly cohesive clusters with high proba- 
bility according to the time-dependent expression pattern 
(Fig. 3). We also made a list of the gene profiles of each clus- 
ter from a total of 310 genes showing similar expression in 
two microarrays by simple t-test analysis (Tables 1-5). On 
the basis of these two different gene lists, we could extract 
genes that not only belonged to the 310 genes but also were 
traced in the cluster analysis. This is the usual method used 
in microarray analysis to narrow a list of genes down to 
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more significant genes. This gene list of 5 different clusters 
will be useful for elucidation of the TGF-|3-mediated molec- 
ular pathways in mammalian cells. 

Microarray data are often criticized owing to poor re- 
producibility [8]. To address this problem, we selected 3 in- 
teresting genes, such as Smad4, Gadd45|3, and TGF-pi, 
and compared the microarray data with RT-PCR data. We 
can place some trust in the microarray data because the 
RT-PCR results were similar (Fig. 4). 

To predict the mechanism by which TGF-(31 induces 
growth inhibition, we examined the altered gene ex- 
pression profile obtained after TGF-pi treatment on the 
basis of known TGF-|3 signal transduction pathways: the 
Smad-dependent pathway and the Smad-independent 
pathway [9,10]. We compared our microarray analysis re- 
sults with known or suggested TGF-|3 signal pathways. We 
found some interesting signal molecules in our microarray 
data. The microarray analysis showed increased ex- 
pression of SARA and the common Smad, Smad4. All of 
these genes are signal molecules regarded as part of the 
Smad-dependent TGF-|3 signal transduction pathway 
[10,11]. From these results, we can assume that the growth 
inhibitory signal transduction of TGF-(31 in bladder cancer 
might be mediated by the Smad-dependent pathway (Fig. 
5). Moreover, the signal intensity of these molecules is in- 
tensified at 48 hours. This result correlated with the cel- 
lular response data in 253J cells after TGF-pi treatment. 
253J cells showed marked growth inhibition after TGF-pi 
treatment and this inhibitory response was pronounced at 
48 hours after treatment. We can assume that the pathway 
is more active at 48 hours after treatment. 

Besides Smad-dependent signal transduction, we found 
7 genes associated with the MAPK pathway. Among these, 
we found increased expression of MAPKK1 and MAPKK4. 
These two genes are critical to activating two different 
MAPK pathways, the ERK-1/2 pathway and the JNK1/2/3 
pathway [12,13]. In the MAPK signal pathway, MAPKK1 
is known to activate Erkl and MAPKK4 is known to acti- 
vate JUN [13,14]. Hence, we can assume that the mecha- 
nism of the growth inhibitory activity of TGF-pi in bladder 
cancer cells may be through activation of the Erkl and JNK 
MAPK signal pathways. 

This work had some limitations. First, the gene list of al- 
tered expression was established from only one bladder 
cancer cell line, 253J. By comparing these data with results 
obtained from another appropriate bladder cancer cell line, 
such as T24 or 5637 cells, we might be able to narrow down 
the gene list. Secondly, the data from the microarray 
should be further confirmed by immunohistochemical 
staining or Western blotting. Nevertheless, our data have 
revealed several TGF-pi-inducible molecular pathways 
involving Smad4 and Gadd45p. Further characterization 
of these molecular pathways may elucidate the growth-in- 
hibitory mechanism of TGF-pi in bladder cancer. 



CONCLUSIONS 

A gene list showing significantly altered expression pro- 
files after TGF-pi treatment was developed according to 
five highly cohesive clusters. The data suggested that the 
growth inhibitory effect of TGF-pi in bladder cancer may 
occur through the Smad-dependent pathway, possibly via 
activation of the Erkl and JNK MAPK kinase pathways. 
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